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Abstract 

We study the t-J model with the inclusion of the so called three-site term which 
comes out from the t/U — > expansion of the Hubbard model. We find that this 
singlet pair hopping term has no qualitative effect on the structure of the pure mean 
field phase diagram for nonmagnetic states. In accordance with experimental data 
on high-T c materials and some numerical studies, we also find wide regions of phase 
coexistence whenever the coupling J is greater than a critical value J c . We show 
that J c varies linearly with the temperature T, going to zero at T = 0. 

PACS: 71.27, 74.20. 



In a recent paper we have presented a new finite temperature phase diagram for the 
nonmagnetic phases (dimer and flux) in the 2-dimensional t-J model below half-filling, 
which shows the existence of broad areas of phase coexistance (dimer-flux or flux- uniform) , 
in accordance with experimental and numerical data on the possibility of separation into 
hole-rich and hole-poor regions @, |3]. [| . 

It is well known that one can obtain the t-J model in the strong coupling regime 
(JJ ^> t) of the one-band Hubbard model: 

H = ~t 1212 c l c i* + U 12 n ^ n U (!) 

(ij) v i 

below half-filling. In (|l|) Latin indices label sites, the annihilation operators for 

a fermion on site i and spin cr =|, j, n io . = cj^c^ and (ij) denotes a sum over ordered 
nearest neighbor (n.n.) pairs. 

The strong coupling limit of (|1|) can be studied either by direct perturbation methods 
|| or by canonical transformations methods || (see also [|7|, || for a detailed review). To 
lowest order in t/U, the Hamiltonian that describes the Hubbard model on the subspace 
with no double occupancies (rii < 1) is given by: 

H = T h + H {1) + H (2) (2) 

where 

T h = - Uj (1 - nis) 4 a c ja (l - n j& ) ; a = -a (3) 

H {1) = E J iJ n in A , J tJ = 4^ f (4 ) 

(ij) u 

H {2) = -— E 12 % t ji ( X ~ c ^ c 3° n i° n n c )t Qt(1 - n w ) (5) 

u m ar 

and the symbol (ij) denotes a summation over unordered pairs of n.n. and (ijl) denotes 
a summation where i ^ I, and both i and I are n.n. to j. The spin operators are defined 

as 

^ = aE <$a?apCif> , (6) 
Z a/3 

with a = (ci, (72, cr 3 ) being Pauli matrices. 

The effective Hamiltonian contains a direct hopping term that gives nonvanishing 
contribution only if the electron jumps from site % to an empty n.n. site j. So we expect 
the effective hopping coefficient t to be reduced by a factor of order 8, S being the doping 
fraction ||. On the contrary which gives rise to a spin dynamics described by the 
standard AFM Heisenberg model, has not to be renormalized, since Si ■ Sj is nonzero only 
if both the sites i and j are singly occupied. is proportional to \t\ 2 and describes 

a double hopping process. Thus, one would be tempted to believe that it has to be 
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renormalized by a factor of 5 2 . However, it is not difficult to see that the three-site 
hopping term H^> gives nonzero contribution only for the situation in which the final 
site I is empty and the intermediate site j is singly occupied. Then this term should be 
renormalized by a factor of 5(1 — 5). Similar conclusions can be reached within different 



approaches pOL O, 12 



For low doping the three-site term is smaller by a factor 5 as compared with the AFM 
term and by a factor \t\/U as compared with the direct hopping term. That is why it is 
usually ignored and one is led to the t-J model fl3|| . However, the Hubbard Hamiltonian 
is often studied in a range of the parameter t/U not so small, especially in applications 
to high T c materials for which one assumes t/U rs 0.1. In this regime the role of the 
three-site term cannot be neglected a priori. 

Numerical works have recently shown that the the three-site term is essential to avoid 
discrepancies between the one-particle spectra of the Hubbard and t-J model in the inter- 
mediate U regime |T^]. According to [14[], only after the inclusion of the three-site term, 



the t-J model can correctly describe the photoemission experiments on cuprate oxides. 

The three-site term describes the hopping of a singlet pair, which can be either rigid or 
with a flip of the spins. Both these processes preserve the local singlet bond and may be 
relevant for superconducting correlations, favoring superconducting states at the expense 
of an antiferromagnetic phase (]lOf| . 

Therefore, according to the discussion above, we approximate the effect of the 
Gutzwiller projector by renormalizing each term of the Hamiltonian (2) with suitable 
factors depending on 5. Our new starting point is hence the Hamiltonian: 

H = ts mi c ^ c i« "tEE c L Cja c\ p c k/3 (7) 

(ij) ot ^ ijkl a/3 

where we have introduced the four fermion matrix 

1 if z n.n. j and i = k, j = I 
■AjiM = \ 5(1 — 5) if i n.n. j, k n.n. /, j = I and i ^ k . (8) 
otherwise 

In JTI], we examined the so-called nonmagnetic phases for the t-J model, first proposed 



by Affleck and Marston |To|, in the mean field approximation and we have presented a 
finite temperature phase diagram (T versus the doping 5). In this paper we shall study 
how the phase diagram changes if we include the three-site hopping term. 

To study the Hamiltonian (|7|) we can parallel the discussion in [|IJ, where we have 
chosen to represent the partition function in the grand canonical ensemble by means of a 
functional integral over Grassmann fields and decoupled the four-fermion interaction by 
means of a Hubbard-Stratonovich transformations over the auxiliary fields Uij(r) = Uj^r). 
In the static approximation Uij(r) = Uij, the partition function reduces to: 
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Z = J [VU*pUij\ exp | K Aa,u Uki \ exp {S eff } (9) 

S eff = -2Tr{log[-G- 1 + /3U]} , (10) 



where 



[G 11' = (iw„ + ^)Mnn' (11) 

[U]L = « , (12) 

[ otherwise 

/x being the chemical potential. 

We consider, as in ||15|| , a pattern of link variables Uij that admits a symmetry for 



translations along the diagonal so that the matrix U depends only on the four parameters 
Uj (j = 1, ■ ■ • , 4) shown in Fig. 1 of As for the pure t-J model, the matrix U can be 
easily diagonalized in momentum space: if we include the three-site term it has eigenvalues 



E k = ± |A fc | (13) 
A fe = xie lkxa + X* 2 e~ ikva + X3e~ ikxa + xle ikya 



with the definitions 



xi = ta + Wi + <$(i-<$)(w 2 * + w 3 + w;) (i4) 

X 2 = t8 + U 2 + S(l-S){Ut+Ul+U 4 ) (15) 
X s = ^ + W 3 + 5(1-5)(^ 1 +ZY*+ZY 4 *) (16) 
X 4 = t6 + U4 + 6{l-S)(u;+U 2 +i4) (17) 

We notice that both the expression for the partition function @ and for the eigenvalues 
of U (0) are formally identical to the one for the pure t-J model we gave in the only 



differences being in the definition of the matrix A and of the parameters x/ s that appear 
in Afe. Thus, to study nonmagnetic phases in the mean field approximation, we can borrow 
from the discussion of [|l| . 

We first discuss in detail the numerical solutions to the self-consistency equations for 
the case J/t = 1. For high temperature and/or high doping the free energy is minimized 
by the uniform phase, characterized by Ui = U G IR + , i = 1, ... ,4. For temperature 
lower than T = J/4 we find two phases: dimer and flux. The former corresponds to two 
solutions of the form IA\ ^ U 2 = U3 = U4 , Uj G (D with \Ui \ ^> \U 2 \ and is stable in the 
lower range of dopings. In the flux phase the link variables Ui are all equal to a complex 
value Ue 1 ^ with (f> 7^ 0. It is the lowest energy state in a range of doping that roughly 
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Figure 1: 

Mean field phase diagram for the t- J model (open circles) 
and with the inclusion of the three-site term (solid triangles) . 



corresponds to that of high-T c superconductivity. The pure mean field diagram is shown 
in Fig. 1 (solid triangles). 



Since the hole mobility is enhanced by the addition of the three-site hopping term, 
we expect disordered phases be favored with respect to the ordered ones. Moreover, the 
three-site term being proportional to 5(1 — 5), this effect should increase with doping 
in the region we considered. This is in fact what we find in our mean field calculation. 
Fig. 1 shows the shifts of the phase boundaries for the usual t-J model (open circles) 
after the inclusion of the three-site term (solid triangles). The region of stability of the 
dimer phase shrinks in favor of flux phase, while the uniform configuration extends to 
considerably lower dopings at the expense of the flux one. 

As we have observed also in for T = our mean field calculations are consistent 
with the ones carried out within a slave boson formalism in [1^, I7|. In particular, we 



agree with [17[ in finding a first order phase transition between the dimer and the flux 



phases. Contrary to ||17|| , we find that the flux-uniform phase transition is of the second 
order for any temperature. This can be seen from either Fig. 3 or Fig. 5 of [|TJ . In Fig. 
3 of IJ it is possible to observe that, increasing S, the angle 0, which distinguishes the 
flux phase from the uniform, decreases continuously from 7r/4 to 0. In Fig. 5 of it is 
easy to see that, at the transition point, there is no jump in the chemical potential, but 
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only a cusp. At this regard, we can also observe that the slope of these curves becomes 
steeper as the temperature decreases. For T = this behaviour gives rise to numerical 
problems and probably this is the reason why in [|17] the authors found a jump in the 



chemical potential and hence concluded that the flux-uniform boundary line describes a 
first order phase transition. 

The diagram of Fig. 1 is constructed without taking into account the possibility of 
coexistence of different phases. Indeed we find two ranges, belonging respectively to the 
dimer and the flux phase, where |f > P|, which clearly indicate an instability to- 
wards phase separation. By means of the Maxwell construction on the chemical potential 
[ffl, we obtain two different regions corresponding to dimer-flux and flux-uniform phase 
coexistence. In the former, the hole-poor (hole-rich) part is in the dimer (flux) phase 
while in the latter the hole-poor (hole-rich) part corresponds to the flux (uniform) phase. 
The final phase diagram is shown in Fig. 2. 




Figure 2: 

The phase diagram after taking into account phase separation. 



Phase separation for nonmagnetic phases is also discussed in |I7 |. In the range of 



dopings relevant to high-T c superconductivity, the authors of |T7j found a wide region of 
phase separation between a half-filling dimer phase and a hole-rich uniform phase, while 
we find a smaller region of phase separation between a half-filling dimer phase and a flux 
phase corresponding to 5 ~ 0.09. This discrepancy is due to a different behaviour of the 
chemical potential as a function of 6. Indeed, within the large- N slave-boson formalism the 
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electronic chemical potential gets shifted by the mean value of the Lagrange multiplier 
that enforces the constraint, while in the approach we adopted such a term is absent. 
Phase separation is a relevant issue because there is a considerable experimental evidence 
that many high-T c superconductors have a regime in which phase separation between an 



AFM hole-poor region and a superconducing hole-rich region occurs [19 . Also, according 



to [^DJ an instability toward phase separation might explain the unusual properties of the 
normal state and superconductivity in copper oxides. This is why the above mentioned 
discrepancies between different approaches deserve a deeper investigation, which we plan 
to describe in a future work. 




Figure 3: 

The region of antiferromagnetic couplings where phase separation occurs. 
The dotted line refers to a calculation without the three-site term. 



We want finally to give our contribution to the open question about the possibility 
that phase separation disappears once J/t is lowered below a critical value J c /t. The 
determination of this quantity is important in order to understand if the t-J model can 
exhibit phase separation in the regime of high-T c superconductors where 0.3 < J/t < 
0.5. Theoretical studies, for T = 0, based on different techniques, lead to contradictory 
results. Using a high temperature expansion the authors of [|l^] find the model phase 
separates only for J/t > 1, i.e. outside the region interesting for superconductivity. 
Similar results have been obtained in some numerical works ||21|| . On the contrary other 
authors find that phase separation occurs for all values of J/t. In particular this has been 
done by means of quantum Monte Carlo calculations [p2| , a combination of numerical 
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and analytical results |2|, |3| and the above mentioned slave-boson approach ||, [T7[]. We 
calculated numerically the critical value J c /t for several temperatures (T = included) 
by observing the disappearance of the instability regions where || > 0. 

The results, plotted in Fig. 3, show that J c /t varies linearly with the temperature and 
vanishes for T = 0, both for the standard t-J model (dotted line) and with the inclusion 
of the three-site term (solid line). Thus, at zero temperature we find phase separation 

22fl . In addition, we 
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for all values of the interaction strength, in accordance with 
can also conclude that at finite temperature there is always a nonzero critical value J c /t 
below which phase separation disappears. In particular there is no phase separation in 
the region interesting for superconductivity for temperatures higher than T « 0.07t. 

Lastly, we observe that in one dimension some exact diagonalization and Monte Carlo 
results 0| suggest that the inclusion of the three-site term suppresses phase separation, 
which is instead present when such a term is not added |24|, [25|] . This seems not to be the 
case in two dimensions, at least at mean field level. Indeed, our calculations show that 
phase separation is still present and the effect of the three-site term is not so dramatic. 



Acknowledgments 

We are indebted to G. Morandi for many valuable suggestions. We are also grateful 
to D. Baeriswyl, F. Napoli and G. C. Strinati for helpful discussions. Finally, we would 
like to thank the referee for suggesting us to look for the occurrence of phase separation 
also at small values of J/t. 



8 



References 



[1] E. Ercolessi, P. Pieri and M. Roncaglia, Phys. Lett. A225 (1997) 331-340 

[2] E. Dagotto and J. Riera, Phys. Rev. Lett. 70 (1993) 682; 
A. Angelucci and S. Sorella, Phys. Rev. B47 (1993) 8858; 
M. Marder, N. Papanicolaou and G.C. Psaltakis, Phys. Rev. B41 (1990) 6920; 
A. Auerbach, B. E. Larson, Phys. Rev. B 43 (1991) 7800; 

[3] V.J. Emery, S.A. Kivelson and H.Q. Lin, Phys. Rev. Lett. 64 (1990) 475; 
S.A. Kivelson, V.J. Emery and H.Q. Lin, Phys. Rev. B42 (1990) 6523 

[4] M. Grilli, R. Raimondi, C. Castellani, C. Di Castro, G. Kotliar, Phys. Rev. Lett. 67 
(1991) 259 

[5] V.J. Emery, Phys. Rev. B14 (1976) 2989; 

[6] M. Aquarone, "Hubbard Correlations in a Single Band", in Physics of Metals, eds. 
E.S. Giuliano and C. Rizzuto, (World Scientific, 1988) and references therein; 
K.A. Chao, J. Spalek and A.M. Oles, J. Phys. CIO (1977) L271; 
K.A. Chao, J. Spalek and A.M. Oles, Phys. Rev. B18 (1978) 3453; 
C. Gros, R. Joynt and T.M. Rice, Phys. Rev. B36 (1987) 381; 
H. Primas, Helv. Phys. Acta 34 (1961) 331; 

H. Primas, Rev. Mod. Phys. 35 (1963) 710 

[7] A. P. Balachandran, E. Ercolessi, G. Morandi and A.M. Srivastava, Hubbard Model 
and Anyon Superconductivity, World Scientific (1990) 

[8] E. Galleani d'Agliano, G. Morandi and F. Napoli, in High Temperature Superconduc- 
tivity, M. Acquarone Ed., World Scientific (1996) 

[9] G. Baskaran, Z. Zou and PW. Anderson, Solid State Comm. 63 (1987) 973 

[10] J. Spalek, Phys. Rev. B37 (1988) 533 

[11] G. Kotliar and J. Liu, Phys. Rev. B38 (1988) 5142 

[12] M. Grilli and G. Kotliar, Phys. Rev. Lett. 64 (1990) 1170 

[13] J.E. Hirsch, Phys. Rev. Lett. 54 (1985) 1317; 

C. Gros, R. Joynt and T.M. Rice, Phys. Rev. B36 (1987) 8190 

[14] H. Eskes and R. Eder, Phys. Rev. B54 (1996) R14226 

[15] J.B. Marston and I. Affleck, Phys. Rev. B39 (1989) 11538 ; 

I. Affleck and J.B. Marston, Phys. Rev. B37 (1988) 3774 



9 



[16] M.U. Ubbens and P.A. Lee, Phys. Rev. B46 (1992) 8434; 
X.G. Wen and P.A. Lee, Phys. Rev. Lett. 76 (1996) 503; 
D.N. Sheng, Z.B. Su and L. Yu, Phys. Rev. B42 (1990) 8732 

[17] M. GriUi, C. Castellani and G. Kotliar, Phys. Rev. B45 (1992) 10805 

[18] W. O. Putikka, M. U. Luchini, T. M. Rice, Phys. Rev. Lett. 68 (1992) 538 

[19] M.F. Hundley, J.D. Thompson, S.W. Cheong, Z. Fisk and J.E. Schirber, Phys. Rev. 
B41 (1990) 4062; 

J.H. Cho, F.C. Chou and D.C. Johnston, Phys. Rev. Lett. 70 (1993) 222; 

J.M. Tranquada, B.J. Sternlieb, D.J. Axe, Y. Nakamura and S. Uchida, Nature 375 

(1995) 561 

[20] A. Perali, C. Castellani, C. Di Castro and M. Grilli, Phys. Rev. B54 (1996) 16216 

[21] P. Prelovsek and X. Zotos, Phys. Rev. B47 (1993) 5984; 
D. Poilblanc, Phys. Rev. B52 (1995) 9201 

[22] C. S. Hellberg and E. Manousakis, |cond-mat/9611195 

[23] B. Amnion, M. Troyer and H. Tsunetsugu, Phys. Rev. B52 (1995) 629 

[24] Masao Ogata, M. U. Luchini, S. Sorella and F. F. Assaad, Phys. Rev. Lett. 66 (1991) 
2388 

[25] F. F. Assaad and D. Wiirtz, Phys. Rev. B44 (1991) 2681 



10 



